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Since  October  1,  1991  experimental  and  theoretical  research  has  been  conducted  by 
Prof.  Min-Chang  Lee  and  his  students  at  BU  and  MIT.  This  research  work  is  aimed 
at  investigating  the  ionospheric  plasma  disturbances  which  can  affect  significantly 
the  radio  wave  propagation  In  communications  and  space  surveillance.  The  research 
topics  which  have  been  Investigated  include:  (1)  A  source  mechanism  leading  to  the 
symmetric  lower  hybrid  sidebands  and  a  low-frequency  mode  in  the  upper  atmosphere, 
(2)  Characteristics  of  lightning-induced  plasmas,  (3)  Radio  wave-produced  plasmas 
and  effects  on  radio  communications,  (4)  Plasma  turbulence  and  formation  of  field- 
aligned  density  fluctuations  as  Ionospheric  ducts. 
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Since  October  1,  1991  experimental  and  theoretical  research  has  been  conducted  by 
Prof.  Min-Chang  Lee  and  his  students  at  BU  and  MIT  under  the  sponsorship  of  Lt. 
Col.  James  G.  Stobie’s  office  at  AFOSR  (grant  ^  AFOSR-92-0001  and  F49620-92-J- 
0297).  This  research  work  is  aimed  at  investigating  the  ionospheric  plasma  disturbances 
which  can  affect  sigmficantly  the  radio  wave  propagation  in  communications  and  space 
surveillance.  The  research  topics  which  have  been  investigated  include:  (1)  A  source 
mechanism  leading  to  the  symmetric  lower  hybrid  sidebands  and  a  low-frequency  mode 
in  the  upper  atmosphere,  (2)  Characteristics  of  lightning-induced  plasm2is,  (3)  Radio 
wave-produced  plasmas  and  effects  on  radio  communications,  (4)  Plasma  turbulence  and 
formation  of  field-aligned  density  fluctuations  2is  ionospheric  ducts.  Progress  in  these 
research  subjects  are  briefly  described  below. 

(1)  A  low-frequency  (~  SOORr)  mode  in  conjunction  with  symmetric  sidebands,  dis¬ 
played  by  approximately  500  Hz  off  the  carrier  frequencies  of  11.90  and  12.65  kHz,  has 
been  observed  in  experiments  conducted  with  the  Aureol  3  satellite  at  altitudes  of  1200 
km  above  the  Komsomolsk-on-Amur  Alpha  Station  (50®iV,  135®E)  (Sotnikov  et  al.,  1991 
and  references  therein).  These  observations  are  examined  and  interpreted  with  an  alter¬ 
native  mech2uiism  based  on  the  theory  of  Lee  and  Kuo  [1984].  The  proposed  mechanism 
suggests  that  the  injected  VLF  wave  parametrically  excites  a  purely  growing  mode  and 
both  the  Stokes  and  anti-Stokes  sidebands  of  the  lower  hybrid  waves.  The  Doppler-shifted 
frequencies  of  these  modes  resulting  from  the  satellite  motion  adequately  account  for  both 
the  sidebands  and  the  low-frequency  mode  observed  in  the  Aureol  3  experiments.  This 


1 


^4  fv'O'/  '^'”9 


work  will  be  published  in  an  JGR  paper  (1992)  by  Dalkir,  Lee,  Groves,  and  Kuo. 

(2)  The  scattering  of  radio  waves  from  lightning-induced  pleismas  is  investigated,  the 
plasmas  are  modeled  as  long,  thin,  perfectly  conducting  vertical  channels  with  surface 
density  irregularities  on  a  scale  much  smaller  than  the  total  channel  length,  the  theoreti¬ 
cally  predicted  wavelength  dependence  of  the  radar  reflectivity  is  compared  to  experimen¬ 
tal  results  obtained  with  the  MIT  C-band,  S-band,  and  the  Millstone-Hill  UHF  radars. 
It  is  found  that  the  rough  surface  scattering  model  can  adequately  account  for  the  1/A^ 
wavelength  dependence  of  the  radar  reflectivity  observed  in  the  MIT-Millstone  Hill  ex¬ 
periments,  if  the  density  fluctuations  of  the  lightning-induced  plasmas  have  a  power-law 
type  spectrum.  The  results  of  this  work  has  been  reported  in  a  paper  submitted  recently 
to  JGR  for  publication  by  Dalkir,  Lee,  and  Williams. 

(3)  The  scattering  of  radio  waves  from  perfectly  conducting,  rough  plasma  surfaces 
is  investigated.  Cylindrical  plasma  channnel  created  by  lightning  discharges  and  disk¬ 
shaped  plasmas  generated  by  a  single  radio  beam  are  examined.  The  surface  density 
irregularities  of  the  plasmas  are  modelled  as  random  perturbations  of  the  constant  den¬ 
sity  background,  characterized  by  either  a  Gaussian  or  power-law  type  spectrum.  The 
theoretical  wavelength  dependence  of  the  radar  reflectivity  of  the  lightning  plasmas  is 
compared  to  experimental  results  obtained  with  the  MIT  C-band,  S-band,  and  the  Mill¬ 
stone  UHF  radars.  The  rough  surface  scattering  model  can  adequately  account  for  the 
1/A^  wavelength  dependence  of  the  radar  reflectivity  observed  in  the  MIT-Millstone  Hill 
experiments,  if  the  surface  density  irregularities  have  a  power-law  type  spectnun.  The 
efficiency  of  reflection  of  radio  waves  from  an  artificially  generated  plasma  in  both  mono¬ 
static  and  bistatic  radar  geometries  is  discussed.  In  radio  conununication  and  space 
surveillance  applications,  the  surface  density  irregularities  may  significantly  adfect  the  re¬ 
flectivity  of  the  artificial  pleisma  patch,  depending  on  the  power  spectrum  of  the  plasma 
density  irregularities.  The  results  of  the  work  are  presented  in  a  paper  submitted  to  the 
technical  journal  titled  “Radio  Science”for  publication  by  Dalkir  and  Lee. 

(4)  In  our  recent  ionospheric  heating  experiments  at  Arecibo,  Puerto  Rico,  we  in¬ 
vestigated  the  spectral  characteristics  of  HF  enheinced  Langmuir  waves  among  other 
phenomena.  The  obliquely  propagating  HF  waves  were  transmitted  from  a  tilted  heater 
in  experiments  during  August  5  -  11,  1992,  while  the  vertically  propagating  HF  waves 
were  employed  in  experiments  during  September  17  -  22,  1992.  In  the  “oblique  heating” 
experiments,  the  main  bezun  of  the  HF  heater  waves  was  tilted  by  12  ®  to  the  north, 
while  the  first  sidelobe  pointed  13  ®  to  the  south.  The  radiated  power  of  the  first  side- 
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lobe  was  estimated  to  be  only  slightly  different  (~  \dB)  from  that  of  the  mean  beeun 
(D.B.  Muldrew,  personal  communication,  1992).  Thus,  the  illuminated  ionospheric  re¬ 
gion  was  presumably  much  greater  than  that  with  the  vertically  transmitted  waves  from 
the  heater.  The  cascading  spectra  of  the  PDI  (parametric  decay  instabilities)  produced 
Langmuir  waves  were  observed  in  the  vertical  heating  experiments  as  expected.  The  radar 
measurements  of  the  HF  enhanced  plasma  lines  in  the  oblique  heating  experiments  zJso 
showed  the  excitation  of  PDI;  however,  the  cascading  spectrum  was  not  seen.  Rather,  a 
broad,  highly  structured  spectrum  separated  from  the  PDI  produced  spectrum  by  15  kHz 
or  so  with  higher  power  was  generated.  It  indicates  that  other  process(es)  excited  by  the 
obliquely  propagating  heater  waves  can  compete  with  the  PDI.  Our  comparative  studies 
of  Langmuir  waves  produced  by  the  vertically  and  obliquely  transmitted  HF  heater  are 
aimed  at  invesitigating  the  source  mechanisms  responsible  for  the  generation  of  HF  en- 
haced  Langmuir  waves.  Pareimetric  instabilities  and  nonlinear  scattering  processes  will 
be  examined  as  the  potential  mechanisms.  This  work  will  be  presented  by  Lee,  Sanchez, 
Kuo,  Moriarty,  Sulzer,  and  Groves  at  the  National  Radio  Science  Meeting  of  URSI  to  be 
held  in  Boulder,  Colorado  during  J  anuary  6-9,  1993. 

These  reseeirch  projects  lead  to  the  publications  and  presentation  of  the  following 
papers  and  students’  theses; 

(A)  Publications  in  Refereed  Journals: 

(1)  Y.R.  Dalkir,  M.C.  Lee,  K.M.  Groves,  and  S.P.  Kuo,  “A  mechanism  responsible 
for  the  observation  of  symmetric  lower  hybrid  sidebeinds  and  a  low  frequency  mode  in 
the  upper  ionosphere”,  accepted  for  publication  in  the  Journal  of  Geophysical  Research, 
1992. 

(2)  Y.R.  Dalkir,  M.C.  Lee,  E.R.  Williams,  “Radio  reflectivity  of  lightning  induced 
plasmas”,  submitted  to  the  Journal  of  Geophysical  Research  for  publicationh,  1992. 

(3)  Y.R.  Dalkir,  and  M.C.  Lee,  “Radar  studies  of  lightning  induced  plasmas  with 
potential  applications  to  radio  communications  and  space  surveillance”,  submitted  to 
Radio  Science  for  publication,  1992, 

(B)  Graduate  Students’  Theses: 

K.D.  Vilece,  “ExperimentaJ  study  of  anti-Stokes  Langmuir  modes  and  particle  pre¬ 
cipitation  produced  by  RF  heating  of  ionospheric  plasmas”.  Master  Thesis,  Department 
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of  Nuclear  Engineering,  MIT,  1992. 

L. A.  Sanchez,  “Studies  of  artificial  spread  F  phenomenon  with  tilted  HF  heater  at 
Arecibo,  Puerto  Rico”,  Master  Thesis,  Department  of  Electrical,  Computer,  and  Systems 
Engineering,  Boston  University,  in  progress,  1992. 

(C)  Presentations  in  Conferences: 

(1)  “Nonlinear  interactions  of  VLF  waves  with  ionospheric  plasmas”,  by  M.C.  Lee,  K.M. 
Groves,  Y.R.  Daikir,  and  S.P.  Kuo,  presented  at  the  URSI-  National  Radio  Science 
Meeting,  University  of  Colorado,  Boulder,  Colorado,  January  7-10,  1992. 

(2)  “Anti-Stokes  Langmuir  modes  produced  by  high  power  HF  waves  at  Arecibo,  Puerto 
Rico”,  by  M.C.  Lee,  K.D.  Vilece,  K.M.  Groves,  and  S.P.  Kuo,  presented  at  the  URSI- 
National  Radio  Science  Meeting,  University  of  Colorado,  Boulder,  Colorado,  January 
7-10,  1992. 

(3)  “Frequency-upshifted  Langmuir  modes  produced  by  RF  heating  of  ionospheric  plas¬ 
mas”,  by  M.C.  Lee,  K.D.  Vilece,  Y.R.  Daikir,  K.M.  Groves,  and  S.P.  Kuo,  presented 
at  the  1992  IEEE  International  Conference  on  Plasma  Science,  Tampa,  Florida,  Jime 
1-3,  1992. 

(4)  “Experimental  studies  of  nonlinear  wave  propagation  Eutid  interactions  with  magne¬ 
tized  plasmas”,  by  M.C.  Lee,  D.T.  Moriarty,  M.  Onozuka,  K.D.  Vilece,  R.R.  Parker, 
and  S.P.  Kuo,  presented  at  the  1992  IEEE  Intemationed  Conference  on  Plasma  Sci¬ 
ence,  Tampa,  Florida,  June  1-3,  1992. 

(5)  “Satellite  measurements  of  symmetric  lower  hybrid  sidebands  and  a  low  frequency 
mode  during  VLF  wave  injection”,  Y.R.  Daikir,  M.C.  Lee,  K.M.  Groves,  and  S.P. 
Kuo,  presented  at  the  1992  International  Beacon  Satellite  Symposium,  Massachusetts 
Institute  of  Technology,  Czunbridge,  Massachusetts,  July  6-10,  1992. 

(6)  “Generation  of  artificial  ionospheric  ducts  for  guided  VLF  wave  propagation”,  by 

M. C.  Lee,  S.  Beisu,  K.M  Groves,  K.L  Koh,  R.  Singh,  M.P.  Sulzer,  and  K.D.  Vilece, 
presented  at  the  1992  Intemationiil  Beacon  Satellite  Symposium,  Massachusetts  In¬ 
stitute  of  Technology,  Cambrdige,  Massachusetts,  July  6-10,  1992. 

(7)  “Laboratory  simulation  of  anomalous  absorption  of  radio  waves  in  turbulent  magne¬ 
toplasmas”,  by  D.T.  Moriarty,  M.C.  Lee,  M.  Onozuka,  R.R.  Parker,  and  K.D.  Vilece, 
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presented  at  the  1992  International  Beacon  Satellite  Symposium,  Massachusetts  In¬ 
stitute  of  Technology,  Cambrdige,  Massachusetts,  July  6-10,  1992. 

(8)  “Analysis  of  wave  scattering  off  a  two  dimensional  plasma  density  irregularities  in 
the  ionosphere”,  by  A.Y.  Ho,  S.P.  Kuo,  and  M.C.  Lee,  presented  at  the  1992  Inter¬ 
national  Beacon  Satellite  Symposium,  Massachusetts  Institute  of  Technology,  Cam¬ 
brdige,  Massachusetts,  July  6-10,  1992. 

(9)  “Scintillation  effects  on  incoherent  scatter  radar  observations  of  enhanced  ion  acoustic 
waves”,  by  K.M.  Groves,  M.C.  Lee,  J.C.  Foster,  and  F.T.  Djuth,  presented  at  the  1992 
International  Beacon  Satellite  S3miposium,  Massachusetts  Institute  of  Technology, 
Cambrdige,  Massachusetts,  July  6-10,  1992. 
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A  Mechanism  Responsible  for  the  Observation  of  Symmetric  Lower  Hybrid 
Sidebands  and  a  Low  Frequency  Mode  in  the  Upper  Ionosphere 


Y.  R.  Dalkir  and  M.  C.  Lee 

Plasma  Fusion  Center,  Massaekuseits  Institute  of  Technology,  Cambridge,  Massachusetts 


K.  M.  Groves 


Ionospheric  Physics  Division,  Phillips  Laboratory,  Hanscom  Air  Force  Base,  Massachusetts 


S.  P.  Kuo 

Weber  Research  Institute,  Polytechnic  University,  Farmingdale,  New  York 


A  low-frequency  (~500  Hz)  mode  in  conjunction  with  symmetric  sideb&nds,  displaced  by  ap¬ 
proximately  500  Hz  off  the  carrier  frequencies  of  11.90  and  12.65  kHz,  has  been  observed  in 
experiments  conducted  with  the  Aureol  3  satellite  at  altitudes  of  1200  km  above  the  Komsomolsk- 
on-Amur  Alpha  Station  (50*N,  135‘’E)  (Sotnikov  et  al.,  1991  and  references  therein).  These 
observations  are  examined  and  interpreted  with  an  alternative  mechanism  based  on  the  theory  of 
Lee  and  Kuo  [1984].  The  proposed  mechanism  suggests  that  the  injected  VLF  wave  parametri¬ 
cally  excites  a  purely  growing  mode  and  both  the  Stokes  and  anti-Stokes  sidebands  of  the  lower 
hybrid  waves.  The  Doppler-  shifted  frequencies  of  these  modes  resulting  from  the  satellite  motion 
adequately  account  for  both  the  sidebands  and  the  low-frequency  mode  observed  in  the  Aureol  3 
experiments. 


1.  Introduction 

Symmetric  electrostatic  sidebands  displaced  by  approxi¬ 
mately  500  Hz  above  and  below  the  carrier  VLF  frequency 
have  been  observed  by  Tanaka  et  al,,  [1987,  and  refer¬ 
ences  therein]  and  others  in  experiments  conducted  at  the 
Komsomobk-on-Amur  Alpha  station  (50®N,  135®E).  The 
sidebands  were  recorded  by  the  Aureol  3  satellite  at  L  - 
values  between  1.7  and  2.0,  at  an  altitude  of  1300  km  above 
the  transmitter,  which  was  operated  at  11.90  and  13.65  KHz. 
The  sidebands  were  observed  only  in  conjunction  with  a  low- 
frequency  mode,  located  at  approximately  500  Hz.  In  an  in¬ 
teresting  follow-up  paper,  Sotnikov  et  al,  [1991]  suggest  the 
low-frequency  mode  b  naturally-occurring  ELF  hiss  simi¬ 
lar  to  that  observed  on  other  low-altitude  satellites  [Gurnett 
and  Bums,  1968;  Afosier,  1971;  Muzzio  and  Angerami,  1972] 
which  couples  nonlinearly  carrier  wave  to  produce  the  ob¬ 
served  sidebands. 

The  purpose  of  thb  paper  b  to  dbcuss  a  different  mech- 
anbm  which  can  also  account  for  the  prominent  features 
of  the  sidebands  and  the  low-  frequency  modes  reported 
by  Tanaka  et  al.,  [1987,  and  references  therein].  In  thb 
proposed  mechanbm,  based  on  the  theory  of  Lee  and  Kuo 
[1984],  the  low-frequency  mode  b  not  a  naturally  occur¬ 
ring  electromagnetic  (em)  emission,  but  b  a  field-aligned 
purely  growing  mode  excited  parametrically  by  the  injected 
VLF  wave.  The  symmetric  sidebands  are  the  Stokes  and 
anti-Stokes  components  of  the  excited  lower  hybrid  waves 
(LHW).  The  theory  shows  that  the  frequency  shifts  of  the 
sidebands  are  comparable  in  magnitude  to  the  observed  fre- 


quency  of  the  low  frequency  mode,  which  is  in  good  agree¬ 
ment  with  the  results  reported  by  Sotnikov  et  al.  [1991]  and 
Tanaka  et  al.  [1987]. 

The  essential  difference  between  the  mechanism  proposed 
in  this  paper  and  that  outlined  by  Sotnikov  et  al.  [1991] 
is  the  interpretation  of  the  low-frequency  mode:  Sotnikov 
et  al.  [1991]  suggest  that  it  is  a  real  frequency,  small  wave 
number  mode;  by  contrast,  we  propose  that  it  has  a  large 
wave  number,  and  that  the  observed  frequency  is  a  Doppler 
effect  arising  from  the  satellite  motion.  The  electric  and 
magnetic  fields  of  the  emissions  were  measured  by  means  of 
spheres  located  on  insulated  booms  on  the  Aureol  3  space¬ 
craft;  this  equipment  cannot  differentiate  a  real  frequency 
from  a  Doppler  shift  effect.  Thus  the  possibility  of  a  field- 
aligned  purely  growing  density  striation  giving  rise  to  the 
observed  low-frequency  mode  cannot  be  ruled  out  [Sotnikov 
et  al.,  1991].  Our  proposed  mechanism  explains  not  only 
the  generation  of  the  LHW  sidebands  but  predicts  the  con¬ 
comitant  occurrence  of  this  low-frequency  mode. 

Section  2  describes  the  highlights  of  the  theory,  and  its 
applicability  to  Tanaka  et  al.  ’t  [1987]  observations.  The  fre¬ 
quency  shift  of  the  sidebands,  the  threshold  incident  electric 
field  for  the  mechanism,  and  the  growth  rate  of  the  instabil¬ 
ity  are  estimated  in  this  section.  Section  3  discusses  some 
further  characteristics  of  the  low-frequency  mode,  including 
its  electrostatic  or  electromagnetic  character  and  its  wave 
vector  angle,  that  are  relevant  to  the  analysis.  Conclusions 
follow  in  section  4. 


2.  Proposed  Mechanism 

Lee  and  Kuo  [1984]  propose  that  symmetric  lower  hybrid 
sidebands  can  be  parametrically  generated  by  an  intense 
VLF  wave  in  the  ionosphere  and  magnetosphere.  In  addi¬ 
tion,  a  field-aligned  purely  growing  mode  with  a  wave  vector 
comparable  in  magnitude  to  those  of  the  sidebands  is  excited 
concomitantly  at  the  expense  of  the  VLF  pump  wave.  The 
frequencies  of  these  waves  as  observed  by  the  satellite  arise 
from  the  Doppler  shift:  /  =  k  •  V«/2t  where  V,,  the  satel¬ 
lite  velocity,  is  approximately  8  km/s,  and  k  represents  the 
wave  vectors  of  the  excited  wave  modes. 

VLF  waves  injected  into  the  ionosphere  from  a  ground 
station  change  to  circularly  polarized  ducted  whistler  waves 
with  wave  vector  ko  =  koOs  and  electric  field  E  =  (it  -f- 
i^)  £0  exp[>'(koz  —  Mo<)],  where  uo  is  the  wave  frequency, 
and  the  z  axis  of  a  Cartesian  coordinate  system  has  been 
taken  along  the  geomagnetic  field.  This  pump  wave  may 
then  excite  the  Stokes  and  anti-Stokes  sidebands  of  lower 
hybrid  waves  [(k_,  u;-)  and  (k^.,u).t.),  respectively]  together 
with  a  purely  growing  or  zero-frequency  mode  (k«,  w,)  ac¬ 
cording  to  the  following  wave  vector  and  frequency  matching 
relations: 


uf+  —  w#  =  u^o  —  4"  (o. 


(1) 


(2) 


k+  —  eik,  =  e,ko  =  k_  +  Cxfc. 

Under  the  assumption  that  the  low-frequency  mode  is 
purely  growing  (w,  =  t-y)  and  that  its  wave  vector  is  per¬ 
pendicular  to  the  geomagnetic  held  (k<  =  ijUx),  equations 
(1)  and  (2)  reduce  to 

=  wo  +  if 


k^  —  Gxko  i  Gxkf 

Furthermore,  k,  »  ko,  and  thus  k±  k,. 

The  magnitude  of  the  wave  number  of  the  sidebands  can 
be  calculated  by  substituting  the  parameters  relevant  for  the 
Aureol  3  experiments  as  provided  by  Sotnikov  et  al.  [1991] 
into  the  lower  hybrid  wave  dispersion  relation: 

Re(w±)  =  Wo  ~  wpi(l  -I- 

where  Wpe/fie  is  the  ratio  of  the  electron  plasma  frequency 
to  the  electron  gyrofrequency.  At  1200  km,  Qe  =  2t  • 
850  kHs  for  Bo  =  3  ■  10~*  T,  and  m/A/  is  the  ratio  of  the 
electron  mass  to  the  ion  mass.  The  magnitude  of  the  wave 
vector  ko  is  obtained  from  the  ducted  whistler  wave  disper¬ 
sion  relation  ko  =  (W|m/c)(wo/0«)’^^.  The  resulting  max¬ 
imum  Doppler  shifted  frequency  of  the  sidebands  is  given 
approximately  by  (wo  ±  |ifc,V',|)/2»,  where  wo/2x  is  the  car¬ 
rier  frequency.  The  Doppler  shifted  frequency  is  reduced 
from  this  value  when  the  satellite  moves  across  the  Earth’s 
magnetic  field  at  an  angle  that  is  not  exactly  90*. 

Figure  1  displays  the  results  of  calculations  of  the  max¬ 
imum  sideband  frequency  displacement  over  a  range  of  ion 
concentration  (wp«/fl«  =  0.5  —  0.85)  for  the  two  carrier  fre¬ 
quencies,  11.90  kHz  and  12.05  kHz.  For  an  ion  composi¬ 
tion  that  is  predominantly  H'*',  the  graph  shows  that  large- 
frequency  displacements  are  possible  near  those  ion  concen¬ 
trations  where  the  carrier  frequency  approaches  the  lower 
hybrid  resonance  frequency.  This  occurs  at  approximately 
Wpx/Oe  =  0.74  and  0.83  for  wo/2x  =  11.90  and  12.65  kHz 
respectively.  A  500-Hz  sideband  frequency  displacement  is 
expected  in  the  range  Wpe/n«  =  0.73  —  0.80.  By  contrast,  if 
is  the  dominant  species,  only  small  frequency  shifts  (< 
100  Hz)  arc  possible.  The  actual  ion  composition  of  the  rel¬ 
evant  region  of  the  ionosphere  can  be  estimated  from  data 
compiled  from  experiments  performed  with  the  Explorer  31 
satellite  [Hoffman  and  Dodson,  1980].  The  altitude  of  1200 
km  at  an  L  -value  of  2.0  is  actually  within  the  transition  re¬ 
gion  between  O'*'  to  H*^  as  the  dominant  ion  species.  Where 
H'^  predominates,  the  daytime  ion  concentration  lies  in  the 
range  Wpe/De  ~  0.70  —  0.80.  Referring  again  to  Figure  1,  it 
can  be  seen  that  this  is  the  range  in  which  a  500-Hz  sideband 
displacement  is  possible. 

Furthermore,  the  requirement  that  as  k,  implies  that 
the  purely  growing  mode  should  be  observed  at  a  Doppler 
frequency  of  /,  =  k,  •  V,/2t  as  k±  •  V,/2x  which  is  approxi¬ 
mately  equal  to  the  frequency  shift  of  the  sidebands.  Figure 
2  of  Sotnikov  et  al.  [1991]  indicates  that  the  maximum  oc- 


currence  of  the  lower  sidebands  is  at  500  to  600  Hz  below 
the  transmitter  frequency,  the  maximum  occurrence  of  the 
upper  sidebands  is  500  to  650  Hz  above  the  transmitter  fre¬ 
quency,  and  the  maximum  occurrence  of  the  low  frequency 
modes  is  between  500  and  600  Hz.  Thus  the  data  are  in 
good  agreement  with  the  predictions  of  Lee  and  Kuo  [1984] 
as  illustrated  in  Figure  1. 

Sidebands  will  be  generated  only  if  the  injected  VLF  wave 
is  intense  enough;  Lee  and  Kuo’s  mechanism  requires  an  inci¬ 
dent  electric  field  threshold  of  [Lee  and  Kuo,  1984,  equations 
(21)  and  (22)] 

domain  1;  u)lh[1  -1- (A//Tn)(u*e/c^)(wpj/flj)]‘^*  <  wo  <  Wp, 

E„  =  (V/m) 

€ 

domain  2:  Wpi  <  wo 

Em  =  0.86(fc^t)?,/n,)(m/e)(l -t-(l 

(V/m) 

where  wtH  is  the  lower  hybrid  resonance  frequency,  Wp,  is 
the  ion  plasma  frequency,  is  the  effective  electron  colli¬ 
sion  frequency,  Vte  is  the  electron  thermal  velocity  and  where 
.?  =  [1  -t-  {M/m){ko/k,f]/[l  -  (M/m)(ko/k,f(n,M% 
Parameters  for  the  upper  ionosphere  have  been  estimated  as 
follows:  Vte  —  ~  1.3  •  10^  m/s  for  T«  ~  lOOOK; 

Ue  =  ^59 -t- 4.181og(5j-))l0~®Tr^^*no  ~  10  Hz;  where  no 
has  been  calculated  from  the  electron  plasma  frequency: 
no  =  Wpe«om/e^.  The  relevant  regime  for  an  H+  plasma 
is  domain  1 ,  with  resulting  threshold  electric  fields  of  Em  = 
8  fiV/m  for  wo/2t  =  11.90  kHz  and  between  9  and  10  /iV/m 
for  wo/2jr  =  12.65  kHz  over  the  range  Wp«/n«  =  0.7  —  0.8. 
A  plasma  whose  primary  ion  composition  is  O'*'  falls  in  do¬ 
main  2,  with  Em  =  15  —  20  /iV/m  for  Upe/Qe  =  0.7  —  0.8 
and  wo/2x’  =  11.90  and  12.65  kHz.  The  Komsomolsk-on- 
Amur  alpha  transmitter  was  operated  at  a  power  of  500  kW 
[Tanaka  et  al.,  1987]  which  can  produce  an  incident  electric 
field  of  at  least  a  few  millivolts  pet  meter  at  an  altitude  of 
1200  km  [5ofniitov  et  al.,  1991];  this  clearly  exceeds  the  re- 
q  .  'red  threshold  for  the  sideband  excitation  by  2  to  3  orders 
of  magnitude.  Thus  the  VLF  pump  waves  are  indeed  intense 
enough  to  excite  the  instability  in  the  upper  ionosphere. 

Lastly,  one  can  estimate  the  optimum  growth  rate  of  the 
instability,  given  by  [Lee  and  Kuo,  1984.  equation  (20b)] 

7  ~  l.HVek]vl/Ql)Eo/Em  is-') 

where  Eo  taken  to  be  10  mV/m  is  the  incident  electric  field 
of  the  pump  wave.  For  an  H'*'  plasma  giving  rise  to  a  500-Hz 
shift,  7^2  s'*' ;  by  contrast,  for  an  O'*’  plasma  producing  a 
100-Hz  sideband  displacement,  7  ~  0.03  s“*.  In  the  Aureol 
3  experiments,  the  alpha  transmitter  pulse  duration  was  0.4 
s  [Tanaka  et  al.,  1987],  which  is  sufficient  for  the  500-Hz, 
H'*'-induced  sidebands  to  be  excited,  but  is  too  short  for  the 
100-Hz  ,  O'*  modes  to  be  induced. 


3.  Discussion 


The  calculations  presented  in  section  3  show  that  the  the¬ 
ory  of  Lee  and  Kuo  [1984]  can  successfully  account  for  the 
500-Hz  frequency  shift  of  the  symmetric  sidebands,  as  well 
the  concurrent  low-frequency  mode  observed  at  500  Hz  by 
Tanaka  et  al.  [1987].  Thb  section  considers  a  few  additional 
aspects  of  the  data  that  further  test  the  applicability  of  the 
proposed  mechanism,  i.e.  the  magnetic  characteristics  of 
the  low-frequency  mode,  the  angle  of  its  wave  vector  to  the 
geomagnetic  held,  power  plots  of  the  phenomenon,  and  the 
ion  concentration  of  the  relevant  region  of  the  ionosphere. 

Sotnikov  et  al.  [1991]  suggest  that  the  low  frequency 
modes  at  500  Hz  are  clearly  electromagnetic  (namely,  ELF 
waves),  whereas  the  symmetric  sidebands  are  electrostatic. 
The  data  plots  of  the  magnetic  held  and  electric  held  in¬ 
tensity  versus  frequency  presented  in  Figure  1  of  their  pa¬ 
per,  however,  do  not  seem  to  support  this  suggestion.  (We, 
however,  note  that  other  data  sets  shown  in  Tanaka  et  al. 
[1987]  do  support  the  suggestion  of  Sotnikov  et  al.  [1991].) 
In  the  original  analysis  of  the  data  by  Tanaka  et  al.  [1987], 
the  minimum  electric  (E)  and  magnetic  (B)  held  intensities 
recorded  by  their  equipment  are  assumed  to  correspond  to 
electromagnetic  emissions.  The  index  of  refraction  calcu¬ 
lated  from  these  minimum  levels  is  found  to  be  cB/E  =  23. 
Au  index  of  refraction  much  smaller  than  this  reference  value 
(e.g.  cB/E  =  7)  implies  that  the  corresponding  mode  is 
quasi-electrostatic.  The  ELF  data  analyzed  by  Sotnikov  et 
al.  yielded  a  value  of  cB/E  statistically  on  the  order  of  60 
(Y.  Tanaka,  personal  communication,  1992).  The  data  pre¬ 
sented  in  Figure  1  of  Sotnikov  et  al.  [1991],  however,  have 
an  index  of  refraction  of  approximately  6  at  500  Hz,  which 
indicates  that  that  this  is  an  electrostatic  mode.  Thus  low- 
frequency  modes  at  500  Hz  or  so  of  both  an  electrostatic 
(ES)  and  electromagnetic  (EM)  nature  have  been  observed 
in  the  Aureol  3  experiments.  Furthermore,  ELF  hiss  phe¬ 
nomena  in  this  frequency  range  have  been  observed  to  be 
intensified  at  latitudes  above  VLF  transmitters,  indicating 
these  modes  may  not  be  entirely  naturally  occurring,  but 
likely  to  be  generated  by  the  VLF  waves  injected  into  the 
ionosphere  at  these  latitudes  [Parrot,  1990], 

A  further  requirement  of  the  proposed  mechanism  is  that 
the  density  striations  produced  by  the  low-frequency  mode 
be  field  aligned,  i.e.,  the  wave  vector  (kj)  should  be  per¬ 
pendicular  to  the  geomagnetic  field.  Sotnikov  et  al.  [1991] 
estimate  the  wave  vector  angle  to  be  9?  >  65®.  In  addition, 
Tanaka  et  al.  [1987]  state  that  there  is  no  direct  evidence  to 
rule  out  electrostatic  characteristics  of  the  ELF  emissions, 
such  as  proton  cyclotron  waves  propagating  almost  perpen¬ 
dicularly  to  the  geomagnetic  field.  Thus  the  requirement  of 
a  tight  angle  between  the  low  frequency  wave  vector  and  the 
geomagnetic  field  is  consistent  with  the  observations. 

Data  plots  of  the  relative  power  of  the  transmitter,  the 
ELF  emission/low  frequency  mode  and  the  sidebands  for  the 
Aureol  3  data  are  presented  in  Figure  6  of  Lagoutte  et  al. 
[1989].  The  authors  show  convincingly  that  the  power  of  the 


sidebands  increases  simultaneously  with  the  low  frequency 
mode,  while  the  received  power  of  the  transmitted  VLF  wave 
is  decreased.  This  is  consistent  with  the  proposed  instability 
process,  which  generates  both  the  sidebands  and  the  purely 
growing  mode  at  the  expense  of  the  incident  pump  wave. 

Lastly,  the  mechanism  proposed  in  section  2  requires  an 
ion  concentration  of  H"*"  in  the  relevant  region  of  the  iono¬ 
sphere  to  produce  a  500-Hz  frequency  shift.  The  mecha¬ 
nism  cannot  operate  if  O'*'  is  the  dominant  species,  since 
the  growth  time  of  the  instability  in  this  case  is  approx¬ 
imately  two  orders  of  magnitude  larger  than  the  incident 
pube  duration.  The  altitude  of  1200  km  at  an  L  -value  of 
2  is  actually  within  the  transition  region  of  the  ionosphere 
from  O'*'  to  H'*'  as  the  dominant  ion  species  [Hoffman  and 
Dodson,  1980].  Which  species  b  actually  more  abundant 
at  the  altitude  of  the  satellite  depends  on  various  factors, 
such  as  the  time  of  day  and  the  solar  cycle  [Hoffman  and 
Dodson,  1980].  The  proposed  mechanbm  can  only  gener¬ 
ate  the  observed  sidebands  and  low-frequency  mode  when 
H+  is  dominant;  thb  may  explain  why  the  sidebands  were 
observed  only  occasionally  by  Tanaka  et  al.  [1987]. 

4.  Conclusion 

The  proposed  mechanbm,  based  on  the  theory  of  Lee 
and  Kuo  [1984],  successfully  explains  the  occurrence  of  sym¬ 
metric  electrostatic  sidebands  together  with  a  low-frequency 
mode  generated  by  VLF  waves  propagating  through  the  up¬ 
per  ionosphere  [5otnil;ov  et  al.,  1991,  and  references  therein}. 
The  parametric  excitation  mechanbm,  whereby  both  the 
Stokes  and  anti-Stokes  components  of  lower  hybrid  waves 
are  generated  together  with  a  purely  growing  mode,  in  a 
plasma  where  H'*'  b  dominant,  can  account  for  the  mag¬ 
nitude  of  the  frequency  shift  of  the  sidebands  and  predict 
the  occurrence  of  the  low  frequency  mode.  The  magnetic 
properties  of  thb  ELF  mode,  its  possible  large  wave  vector 
angle  with  respect  to  the  geomagnetic  field,  studies  of  the 
relative  power  intensities  of  the  generated  modes  and  the 
transmitted  wave,  and  the  ion  concentration  in  the  region 
are  all  consbtent  with  the  proposed  mechanbm.  Thus  we 
conclude  that  a  Doppler  shift  effect,  in  contrast  to  the  real 
frequency  modes  proposed  by  Sotnikov  et  al.  [1991],  can 
also  be  a  plausible  mechanbm  responsible  for  both  the  side¬ 
band  structures  and  the  low  frequency  modes  observed  in 
the  Aureol  3  experiments. 
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Figure  1.  Calculated  frequency  (frequency  shift),  due  to  the  pro¬ 
posed  Doppler  shift  effect,  of  the  low  frequency  mode  (the  sym¬ 
metric  sideband  modes)  as  a  function  of  the  ratio  of  the  plasma 
frequency  (uipt)  to  the  electron  cyclotron  frequency  (Oe)  when  (l) 
the  dominant  ion  species  is  H'f'  and  the  incident  wave  frequency 
is  11.90  kHz  (denoted  by  the  dashed  line),  (2)  the  dominant  ion 
species  is  H'*'  and  the  incident  wave  frequency  is  12.65  kHz  (de¬ 
noted  by  the  full  line),  (3)  the  dominant  ion  species  is  and  the 
incident  wave  frequency  is  either  11.90  kHz  or  12.65  kHz  (denoted 
by  the  dash-dot  line)-no  significant  difference  is  found  between  the 
results  for  these  two  frequencies.  At  the  altitude  of  1200  km.  the 
measured  u>p«/n«  was  found  to  be  between  0.7  and  0.8  [Hoffman 
and  Dodson,  1980]. 
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Figure  1.  Calculated  frequency  (frequency  shift),  due  to  the  proposed  Doppler  sluft  effect,  of  the 
low  frequency  mode  (the  symmetric  sideband  modes)  as  a  function  of  the  ratio  of  the  plasma 
frequency  (u/p*)  to  the  electron  cyclotron  frequency  (0«)  when  (1)  the  dominant  ion  species  is 
H+  and  the  incident  wave  frequency  is  11.90  kHz  (denoted  by  the  dashed  line),  (2)  the  domintint 
ion  species  is  H"*"  and  the  incident  wave  frequency  is  12.65  kHz  (denoted  by  the  full  line),  (3) 
the  dominant  ion  species  is  0+  and  the  incident  wave  frequency  is  either  11.90  kHz  or  12.65  kHz 
(denoted  by  the  dash-dot  line)-no  significant  difference  is  found  between  the  results  for  these  two 
frequencies.  At  the  altitude  of  1200  km,  the  measured  urpe/fle  was  found  to  be  between  0.7  and 
0.8  [Hoffman  and  Dodson,  1980]. 
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Abstract 


The  scattering  of  radio  waves  from  lightning-induced  plasmas  is  investigated. 
The  plasmas  are  modeled  as  long,  thin,  perfectly  conducting  vertical  channels 
with  surface  density  irregularities  on  a  scale  much  smaller  than  the  total  channel 
length.  The  theoretically  predicted  wavelength  dependence  of  the  radar  reflectivity 
is  compared  to  experimental  results  obtained  with  the  MIT  C-band.  S-band.  and 
the  Millstone  UHF  radeirs.  It  is  found  that  the  rough  surface  scattering  model  can 
adequately  account  for  the  1/A-  wavelength  dependence  of  the  radar  reflectivity 
observed  in  the  MIT- Millstone  Hill  experiments,  if  the  density  fluctuations  of  the 
lightning-induced  plasmas  have  a  power-law  type  spectrum. 
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1  Introduction 


Lightning  is  a  series  of  short,  complex,  pulsed  electric  dischcirges  that  create  highly 
irregulcu-  plasma  strips  in  the  atmosphere.  A  typical  lightning-induced  plasma  chcuinel 
is  kilometers  in  length  and  centimeters  to  meters  in  diameter,  with  a  complex  dendritic, 
tortuous  structure  [Uman.  1969].  Radar  lightning  studies  have  been  conducted  since 
the  19o0’s  using  X-band  to  L'HF  radio  waves,  and  extensive  data  are  available  in  the 
literature.  In  particular,  the  wavelength  dependence  of  the  radar  reflectivity  is  of  interest, 
since  it  is  thought  to  be  related  to  the  surface  density  inhomogeneities  of  the  plasma 
chzinnels  [Lee  et  ah,  1988). 

This  paper  is  concerned  with  a  radar  lightning  study  conducted  with  the  MIT  S* 
band,  C-band,  and  the  Millstone  UHF  radars.  These  experiments  demonstrated  an 
overall  inverse  square  dependence  between  radar  reflectivity  and  wavelength,  though 
the  results  were  highly  variable  at  any  one  wavelength,  and  likely  to  be  affected  by 
precipitation  maisking  [Williams,  et  al.,  1989].  Previous  theories  of  lightning  beickscatter 
range  from  modeling  the  plasma  as  an  overdense  radar  target,  i.e.  a  perfect  conductor 
[Ligda,  1956;  Ma^zur  et  ah.  1985],  to  considering  it  ^ls  composed  of  underdense  plasma 
channels  of  varying  geometries  [Atlas,  1958;  Hewitt,  1957;  Holmes  et  al.,  1980].  None 
of  the  published  theories,  however,  adequately  accounts  for  the  average  inverse  square 
wavelength  dependence  of  the  radar  reflectivity  of  the  lightning-induced  plasmas. 

Williams  et  al,  [1989]  attempt  to  resolve  these  various  conflicting  theories  by  pre¬ 
senting  a  clear  case  for  an  overdense  plasma,  based  on  observations  of  temp>erature  in 
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thundercloud  currents  and  laboratory  discharges,  together  with  theoretical  calculations 
of  the  temperature  dependence  of  the  plasma  frequency.  They  conclude  that  one  may 
safely  model  the  lightning  plasmas  as  perfectly  conducting,  thin  channels  of  complex 
geometry.  The  failure  of  the  existing  theories  to  account  for  the  wavelength  dependence 
of  the  radar  reflectivity  may  be  due  to  the  neglect  of  the  small  scale  irregularities  of  the 
plasma  strips  in  the  models.  The  present  study  attempts  to  correct  this  neglect  by  in¬ 
vestigating  the  problem  of  diffuse  scattering  of  electromagnetic  waves  from  an  irregtilcu*. 
perfectly  conducting  plasma  surface.  The  small-scale  plasma  density  irregularities  are 
aissumed  to  have  a  power-law  rather  than  a  Gaussian  type  spectrum,  as  suggested  by 
wave  propagation  studies  conducted  in  both  laboratories  and  in  the  ionosphere  [Lee  et 
al.,  1989].  The  theoretical  wavelength  dependence  of  the  reflected  power  in  this  model 
is  found  to  be  in  good  agreement  with  the  observations. 

The  paper  is  organized  as  follows:  Section  2  presents  a  summary  of  the  relevant 
results  obtained  in  the  MIT-Millstone  Hill  experiments;  Section  3  outlines  the  proposed 
theory.  The  assumptions  made  in  the  model  of  radio  wave  scattering  off  a  plasma 
strip  with  density  irregularities  (viz.,  rough  surface)  are  further  discussed  in  Section  4. 
Conclusions  follow  in  Section  5. 
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2  Experimental  Results 


The  radax  lightning  experiments  were  carried  out  with  the  S-band  (11  cm)  and  C-band 
(5.4  cm)  radaxs  on  the  MIT  campus,  and  the  steerable  150  foot  UHF  (68  cm)  radax 
dish  operated  by  the  MIT  Haystack  Observatory  at  Millstone  Hill  [Williams  et  al..  1989: 
Dalkir  et  al.,  1992].  The  radaxs  were  pulsed  for  1  to  3  fis  at  repetition  rates  of  between 
200  to  900  Hz.  The  pulse  echoes  were  recorded  over  a  wide  range,  from  a  few  kilometers 
up  to  one  hundred,  over  different  storm  days.  The  radax  beamwidth  was  1.4°  for  the 
MIT  radaxs,  and  1.0°  for  the  Millstone  UHF  radax.  The  S-band  and  C-band  radaxs  were 
operated  with  horizontal  electric  field,  while  the  UHF  beam  was  circulaxly  polarized. 
.\11  experimental  runs  were  conducted  with  the  radax  beams  fixed  (See  Figure  1). 

A  histogram  of  the  observations,  showing  the  large  variability  of  the  volume  reflec¬ 
tivity  at  any  one  wavelength,  may  be  found  in  Figure  11  of  Williams  et  al.  [1989].  The 
exact  causes  of  this  variability  remaiin  unexplained  at  present,  but  is  likely  to  be  due 
to  the  variation  in  size,  shape  and  location  of  the  lightning-induced  pletsmas  over  the 
course  of  the  experiments.  Reflection  from  the  precipitation  accompanying  the  lightning 
storms,  or  “precipitation  masking”,  may  also  effectively  reduce  the  observed  wavelength 
dependence  of  the  backscattered  pnjwer  [Williams  et  al.,  1989].  The  mean  reflectivity  at 
any  one  wavelength,  however,  clearly  displays  a  1/A^  dependence  (see  Figure  2). 
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3  Theory 


3.1  Radio  Reflectivity  of  Lightning-Induced  Plasma  Channels 

This  section  outlines  the  theory  of  rough  surface  scatter,  cis  applied  to  a  lightning- 
induced  plasma  layer  in  the  atmosphere.  The  analysis  follows  the  method  introduced 
by  M.  A.  Isakovich  (1952),  who  considered  the  scattering  of  sound  waves  from  perfectly 
reflecting  surfaces  with  randomly  distributed  irregularities.  His  method  wcis  extended 
by  P.  Beckmann  and  A.  Spizzichino  (1963)  to  include  electromagnetic  wave  scattering 
from  perfectly  conducting  surfaces  of  varying  degrees  and  types  of  roughness.  The 
following  analysis  exp£uids  these  results  by  appljnng  them  to  the  specific  problem  of 
radar  backscatter  from  lightning-induced  plasmas. 

Consider  the  scattering  of  plane  electromagnetic  waves  with  electric  field  £'ie**'i  *‘~*‘**, 
where  ki,  ^  are  the  wave  vector  and  angular  frequency,  respectively,  and  r  is  the  position 
vector  with  respect  to  the  origin  of  a  cylindrical  coordinate  system  with  basis  vectors 
er,e(^,ez.  The  scattering  plasma  surface  is  most  appropriately  modeled  as  a  long,  thin, 
rough  cylinder,  with  its  axis  of  symmetry  aligned  with  the  z-£ixis  (see  Figure  3).  The 
vector  s  =  [5  -I-  <5s(z,  <A)]er  +  extends  from  the  origin  to  surface  of  the  plasma.  The 
unperturbed  plasma  channel  has  radius  S;  first-order,  random  deviations  from  this  are 
described  by  the  function  ^s(z,  (f>).  The  scattered  wave  is  characterized  by  the  pairameters 
E2,k2  and  w. 

In  this  simplified  model,  only  the  small-scale  density  irregularities  along  any  one 
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section  of  an  approximately  vertical  lightning  channel  are  considered.  The  breinch- 
like  structure  of  the  plasma  channels,  or  their  large-scale  deviations  from  the  surface 
5  =  5  (i.e..  any  background  inhomogeneities),  are  ignored.  Furthermore,  since  the 
overdense  plaisma  channels  are  long  and  thin,  it  is  reasonable  to  assume  that  any  density 
perturbations  will  have  azimuthal  symmetry:  henceforward  we  assume  6s(c.o)  = 

The  analysis  proceeds  by  use  of  the  Helmholtz  integral  which  describes  the  total 
electric  field  at  any  point  P  in  terms  of  the  electric  field  on  the  scattering  surface: 

where  /Pf,  Pf  is  the  distance  from  P  to  the  scattering  point.  Ea  is  the  electric 

field  on  the  surface  ^4,  and  n  is  the  normal  vector  to  the  surf6w:e  at  the  scattering  point. 
To  deal  with  plane  waves,  we  consider  only  the  far  field,  i.e.  the  point  P  is  distant  from 
the  scattering  poirii,  such  that  Arj/i'  «  k2Po  —  k2  •  s  and  ^  =»  exp(iA:2i2o  ~  *^2  •  s)/Po 
(see  Figure  4).  Furthermore,  we  ignore  multiple  reflection  from  the  rough  surface,  so 
that,  in  the  steady  state, 

=  (!  +  «)£, e'M  (2) 

(3) 

where  R  is  the  reflection  coefficient  and  n  is  the  local  normal  to  the  surface  at  the  scat¬ 
tering  point.  The  reflection  coefficient  is  in  general  a  function  of  the  angle  of  incidence 
and  polarization  of  the  incident  waves,  and  the  electrical  properties  of  the  reflecting 
material. 
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Equation  (2)  is  a  good  approximation  to  the  surface  field  when  the  radius  of  curvature 
of  the  rough  surface  is  large  with  respect  to  the  wavelength  of  the  incident  w'ave.  i.e.  the 
surface  has  no  sharp  edges.  The  normail  derivative  of  the  surface  electric  field  is  obtained 
by  considering  the  reflected  and  incident  waves  with  respect  to  the  tangent  plane  to  the 
surface  at  the  scattering  point.  Equation  (3)  is  valid  within  the  ray  approximation,  or 
when  edge  (diffraction)  effects  can  be  ignored.  This  is  the  case  for  L  »  A.  where  L 
represents  the  characteristic  linear  size  of  the  scattering  surface,  and  A  is  the  wavelength 
of  the  incident  wave.  typical  lightning-induced  plasma  cylinder  is  kilometers  in  length 
and  may  be  meters  in  diameter,  and  thus  satisfies  the  condition  Z,  >>  A  in  both  the  Cr 
and  Cz  directions  [Uman,  1969]. 

Substituting  Equations  (2)  and  (3)  into  the  Helmholtz  integral  (Equation  (1))  yields 

E{P)  =  / / [«(ki  -  k2)  •  n  -  (ki  +  k2)  •  (4) 

It  is  convenient  to  consider  the  normalized  electric  field  E  =  E l[iE\e'’^^ /AnRo]-  Fur¬ 
thermore,  k2  =  —  ki,  since  monostatic,  diagnostic  radars  were  used  in  the  experiments. 
Equation  (4)  can  thus  be  rewritten  ais: 

E  =  2Rj  Jkne^^^^dA  (5) 

where  k  =  ki;  R,  the  local  reflection  coefficient  at  each  scattering  point,  has  been 
aproximated  to  be  constant  over  the  entire  surface.  This  is  appropriate  for  the  lightning- 
induced  plasmas,  which  are  almost  perfectly  conducting,  i.e.  \R\  as  1  [Williams  et  al., 
1989]. 
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The  surface  element  and  normal  vector  can  be  written  in  terms  of  the  density  per¬ 
turbation  65(c)  as:  _ 

dA  =  Sdodzyjl  -H  (65')2  (6) 

n  =  (-6s'ez  -t-  er)^l  + 

where  65'  =  d(6s)/dz.  Substituting  these  terms  into  Equation  (5)  yields 

E  =  2R  I  Sdo  I  dz[-kM  +  (7) 

Jo  J-h/2 

where  h  is  the  characteristic  height  of  the  scattering  cylinder. 

After  performing  the  integration  over  o.  and  integrating  the  term  — 
by  peirts,  Equation  (7)  becomes 

E  =  4/2e‘^7r5[(^^^)  _  4/?^  sin(ib,/i)]  (8) 

>■  kr  J~h/2  kr  J 

where  for  simplicity  we  have  assumed  that  the  surface  is  fixed  at  the  edges;  i.e.  6s{h/2)  = 

6s{—h/2)  =  0.  The  contribution  of  the  second  term  to  the  scattered  power  is  small 

compared  to  the  first,  as  may  be  verified  by  retaining  it  throughout  the  analysis;  in  the 

interest  of  simplicity,  it  will  be  neglected. 

The  mean  scattered  power  can  be  expressed  as  <  P  >=  (1/2)  <  EE*  >  /rjo  where 
Tjo  is  the  intrinsic  impedance  of  the  propagation  medium;  tjq  =  1207r  ~  377  (fi)  for  free 
space.  From  Equation  (8),  we  may  write 

<  EE*  >=  +  dzi  >  (9) 

K*  J-hl2  J-h/2 

where  we  have  taken  =  1 . 

To  proceed  further,  one  must  adopt  a  specific  model  for  the  roughness  parameter 
6s(c).  Plasma  irrregularities  created  by  lightning  are  unlikely  to  have  a  periodic  or 
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similarly  deterministic  structure;  rather,  noise-like,  random  perturbations  are  expected. 
We  therefore  model  the  e.xtent  of  the  plasma  in  the  r  direction  cis  a  random  variable 
65(2).  which  is  completely  described  by  a  Gaussian  joint  probability  density  function 
/[6s(2i),  65(22)]  with  standard  deviation  a.  and  a  normalized  autocorrelation  function 
^(1-1  —  -2I)  with  characteristic  length  I  and  index  m: 


/[65(2i), 65(22)]  = 


1 


—  C- 


65(21)2 -2CC(2i)6s(22) +65(22)-, 

«p< - 2^^(1-C^) - > 


C(|2i-22|)=exp(-|2i-22r/r) 


(11) 


Physically,  <t  and  I  characterize  the  height  and  size  of  the  irregularities  of  the 
lightning-induced  plasma.  Any  roughness  profile  can  be  approximated  by  varying  the 
parameters  of  the  joint  probability  distribution  and  the  autocorrelation  function:  a  large 
value  of  a  and  a  small  I  imply  highly  localized  peaks  with  large  density  variations;  by 
contrast,  small  a  and  large  I  imply  generally  low,  wide  irregularities  (see  Figure  3). 
In  order  to  realistically  ignore  multiple  reflections  in  the  model,  we  have  «issumed  the 
absence  of  sharp  peaks  in  the  roughness  distribution  on  the  surface  of  the  perfectly  con¬ 
ducting,  lightning-induced  plasma;  this  requirement  is  satisfied  in  the  regime  /  >>  A, 
I  >  cr.  Previous  authors  have  concentrated  on  Gaussian  spectra  (i.e.  m  =  2  in  Equation 
(11))  [Isakovich,  1952;  Beckmann  and  Spizzichino,  1963];  studies  of  transionospheric 
radio  wave  propagation  have  shown,  however,  that  power-law  type  spectra  provide  the 
more  realistic  model  for  plasma  irreguleirities  [Lee  et  al.,  1989,  and  references  therein]. 
We  therefore  take  m  —  \  for  the  function  given  in  Equation  (11). 
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The  definition  of  the  mean  in  Equation  (9)  can  then  be  rewritten  in  terms  of  the 


probability  density  functions  as: 


<  >=  f  dzi  r  =  x2(2Av. -2A', 

*/— oo  */— oo 


(12) 


where  \2  is  the  two-dimensional  characteristic  function  of  the  joint  probability  distribu¬ 
tion  /(zi^zo)-  Similarly,  the  one-dimensional  characteristic  function  is  defined  ais 


Xi2kr)  =<  >=  r  f{z)e'^->^^^dz  (13) 

where  f{z)  =  (v/2Tcr)~^  exp(— 2-/2<r-).  It  is  useful  to  first  calculate  the  variance, 
<^{^)=  <  EE*  >  —  <  E  X  E  >*,  given  by 

(HE)  =  <izje‘“-l''-'>l[x2(2fc,  -2fc)  -  v(2fc,).x'(2lt,)) 

(14) 

For  a  Gaussian  joint  probability  distribution  function,  \;(2A:r)  =  exp{-2k'^cr-)  and  X2(2A:r,  -2it,.) 
6xp[~4A:^<7^(l  —  C)];  thus  Equation  (14)  becomes 

(i(f:)=4ir2S^i£iS!  (15) 

The  integrand  vanishes  when  C  =  0,  or  when  jzi  —  22I  is  large;  thus  we  may  appoximate 
C  =  exp(— |zi  — 22!)//  «  1  — |2i  — ^2|/I  and  replace  the  integration  limits  for  22  by  > 

22,  where  /  is  the  maximum  range  of  (21  —  22I  for  which  the  first  order  approximation  to 
the  exponential  is  valid.  Thus  Equation  (15)  can  be  expressed  as 

(HE)  =  dz,  £'*'  (16) 
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This  is  a  good  approximation  to  the  variance  insofar  as  the  integrand  in  Equation  (15) 
is  localized  within  ±/;  it  can  be  seen  that  this  is  the  case  for  a  very  rough  surface, 
i.e.;  ik~a-  >>  1.  which  is  the  relevant  regime  for  the  lightning  backscatter  experiments 
considered  here  (see  Section  3.1.1).  Integrating  the  expression  in  Equation  (16)  yields 
the  following  final  result: 


=  4hff' 


<  EE*  >=  diE)+  <E><E>* 
-2e 


(17) 


=  l2fcr/ 


Ak^a^/l  -  i2k 


37^1 


sm(2kj)/k,  +  sm-{k,h)/{k;h) 


where  Re  means  the  real  part;  the  quantity  <  E  ><  E  >*  has  been  calculated  directly 
from  Equations  (8)  and  (13).  Neglecting  terms  on  the  order  of  yields; 

,1.1 

The  term  (k^  4-  k^)  =  k^  —  k^  indicates  that  <  P  >=  0  for  A:  =  k^,-,  i.e.  the  backscat- 
tered  power  is  zero  for  greizing  incidence  of  the  radio  waves  on  the  plasma  cylinders,  as 
expected. 


3.1.1  Comparison  to  Experimental  Results 

Figure  3  illustrates  the  scenario  of  a  typical  experimental  run:  the  transmitting  radar 
may  be  located  tens  of  kilometers  from  the  lightning  bolts,  which  extend  more  or  less 
verti  ::ally  through  the  atmosphere  for  an  average  of  ten  kilometers  for  cloud  to  ground 
lightning  [Uman,  1969].  A  long  cylindrical  pleisma  channel  is  generated  around  the 
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lightning  discharge,  with  first  order  random  irregulcirities  extending  orthogonally  to  the 
surface.  Observations  of  the  luminous  core  of  the  plasma  indicate  typical  diameters 
of  1  cm  or  less,  where  the  lightning  arc  temperature  may  exceed  5000° K  for  tens  of 
milliseconds  [Williams  et  ah.  1989].  This  corresponds  to  electron  densities  of  between 
10^^- 10^°  m~^,  indicating  that  the  core  of  the  plasma  is  overdense  to  incident  radio  fre¬ 
quencies  of  up  to  approximately  80  GHz.  An  electron  density  of  one  order  of  magnitude 
less  would  still  be  overdense  to  the  radio  frequencies  used  ‘  the  MIT-Millstone  experi¬ 
ments.  Thus  the  effective  scattering  region  is  expected  to  extend  to  radii  larger  than  the 
luminous  core,  where  a  cooler,  less  dense  plasma  has  been  formed.  Previously  reported 
photographic  observations  and  theoretical  calculations  predict  plasma  channel  radii  on 
the  order  of  meters  [Uman,  1969,  and  references  therein].  Thus  we  expect  the  effective 
overdense  plasma  cylinders  to  be  kilometers  long  and  meters  in  diameter. 

The  height  (cr)  of  the  surface  density  irregularities  of  these  plasma  cylinders  cannot 
be  directly  observed,  since  they  occur  in  the  non-luminous  outer  regions  of  the  plasma. 
Based  on  the  expected  dimensions  of  the  plasma  tubes,  however,  one  may  reasonably 
estimate  the  standard  deviation  a  to  be  on  the  order  of  tens  of  centimeters,  and  the 
outer  scale  length  of  the  irregularity  spectrum  or  autocorrelation  distance  I  to  be  on 
the  order  of  meters.  This  satisfies  the  required  conditions  that  the  surface  density 
irregularities  be  smooth,  i.e.  I  »  X  and  <t  <  1.  Under  these  conditions,  it  can  be  shown 
that  >  4k^cr*fP,  and  thus  the  meam  b£ickscattered  power  in  Equation  (18)  can  be 
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approximated  as 


<  P  >2= 


S-cr^-(k;  +  k^r- 

^  7r2  21  kl 


X  1/A- 


(19) 


The  last  proportionality  is  true  since  the  angle  of  incidence  {$)  of  the  radio  waves  on  the 
scattering  surfaces  is  not  expected  to  vary  significantly  over  the  course  of  the  experiments 
(see  Figure  3). 

The  full  range  of  validity  of  Equation  (19)  must  be  carefully  examined.  The  result 
may  be  applied  only  where  the  approximation  used  in  evaluating  the  integral  in  Equation 
( 15)  is  justified;  this  is  the  case  for  >>  1-  Equation  (19)  is  valid  only  in  the  limit 

k^  »  Ak^cr’^JP.  Furthermore,  multiple  reflections  from  the  scattering  surface  may  be 
ignored  only  insofar  as  the  density  irregularities  are  smooth;  this  condition  is  met  when 
I  »  \.  The  various  requirements  for  which  Equation  (19)  is  valid  may  therefore  be 
summarized  as  follows:  Ak^ff^  »  1;  kl  »  Ak*(T^/P  and  A  <<  /.  These  bounds  may 
be  written  in  terms  of  the  angle  of  incidence  of  the  radio  waves  cis  4;rcos^  >  A/<t  > 
[cr !  1)  Air  {cos  6 !  tan0);  A/cr  <<  //cr,  where  for  simplicity  we  have  assumed  that  =  0. 
Figures  5a  and  5b  illustrate  these  relations  for  two  different  values  of  (t/ L  It  can  be  seen 
that  the  region  of  validity  of  Equation  { 19)  is  relatively  large  and  sensitive  to  the  exact 
value  of  the  dimensions  of  the  surface  irregularities. 

The  geometry  illustrated  in  Figure  3  indicates  that  large  angles  of  incidence  {6  >  45°) 
are  expected  in  the  experiments.  Intracloud  lightning  discharges  are  far  more  common 
than  cloud-to-ground  lightning;  the  radio  waves  are  more  likely  to  be  scattering  off 
plasma  channels  that  are  generated  high  above  the  ground.  If  we  estimate  B  w  60°,  on 
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average,  Figure  4b  indicates  that  the  analysis  is  valid  only  in  the  regime  0.5  <  A/a  <  6. 
This  condition  is  indeed  satisfied  for  the  estimated  sizes  of  the  density  irregularities  (a  « 
10  cm;  I  ss  2  m). 

The  one-dimensional  model  adopted  in  the  analysis,  i.e.  6s{z,0)  =  <5s(c)  is  appro¬ 
priate  for  the  lightning-created  plasmas,  which  are  long  cylinders  aligned  with  the  z-axis 
(see  Figure  4).  The  random  density  fluctuations  of  the  plasma  in  the  r-direction  are  thus 
expected  to  depend  only  on  the  z-coordinate  of  the  surface,  and  not  vary  significantly 
with  azimuth.  Since  the  length  (h)  oi  the  lightning-induced  plasma  strips  is  very  long 
(on  the  order  of  a  few  kilometers),  the  other  conditions  assumed  in  the  ancdysis  ,  i.e. 
X  «  I  «  h,  are  reasonably  met  for  all  three  frequencies  used  in  the  experiments. 

The  experimental  data  shown  in  Figure  2  are  thus  within  the  regime  of  validity 
of  Equation  (19),  which  predicts  an  inverse  square  fwjwer  law  dependence  between  the 
backscattered  power  and  the  incident  wavelength.  The  actual  line  of  best  fit  for  the 
experimental  data  indeed  show  an  inverse  square  dependence  of  the  power  on  the  radar 
wavelength,  which  is  in  good  agreeement  with  the  theoretical  result  of  Equation  (19). 
It  should  be  stressed  that  the  relation  <  P  >oc  1/A^  in  Equation  (19)  only  holds  for 
pleisma  density  irregularities  with  a  power-law  spectrum  whose  correlation  function  is 
given  by  Equation  (11)  with  m  =  1 .  Both  laboratory  and  ionospheric  studies  of  turbulent 
plaamais  indicate  that  density  perturbations  are  best  modeled  by  these  types  of  spectra 
[Lee,  et  al.,  1989].  A  Gaussian  spectrum  for  the  density  irregularities,  by  contrast,  leads 
to  <  P  >oc  1/ A,  which  does  not  agree  with  the  observations. 
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4  Conclusions 


riu*  vvavelenglh  dependence  of  the  backscattered  power  from  lightning-induced  plasmas 
has  l)een  studied  in  experiments  conducted  with  the  MIT  S-band.  C-band,  and  the 
Millstone  UHF  radars.  The  average  reflected  power  was  observed  to  decrease  with  the 
wavelength  as  1  /  This  behavior  can  be  successfully  explained  by  the  scattering  of  elec¬ 
tromagnetic  waves  from  the  small-scale  density  irregularities  of  the  lightning-induced, 
perfectly  conducting  plasma  channels  in  the  atmosphere.  These  density  irregularities 
are  reasonablv  assumed  to  have  a  power-law  rather  than  a  Gaussian  spectrum.  A  re¬ 
gion  of  validity  for  the  predicted  1/A^  relation  is  believed  to  exist  in  the  experiments, 
w’hich  depends  on  the  expected  size  of  the  density  irregularities  (<r//),  and  the  angle  of 
incidence  of  the  radio  waves  on  the  plasma  channels. 
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Figure  Captions 


Fig.  1.  Schematic  representation  of  lightning  backscatter  experiments.  The  plasma 
cylinders  created  through  air  breakdown  <iround  the  lighning  discharges  is  diagnosed 
with  the  MIT  C-band.  S-band  and  Millstone  UHF  radars.  For  simplicity,  only  one  light¬ 
ning  discheirge  is  shown:  in  actuality,  the  experiments  were  carried  out  separately  in 
different  electrical  storms. 


Fig.  2.  Summary  of  experimental  results.  The  average  volume  reflectivity  at  A  =  5.4, 
1 1  and  68  cm  is  plotted.  The  numbers  in  parentheses  indicate  the  number  of  data  points 
collected.  The  wavelength  dependence  pte<Ucted  by  the  theory  of  long  thin  conductors, 
and  underdense  plasma  channels  is  also  indicated  on  the  figure.  Neither  existing  model 
fits  the  data  of  the  present  study,  which  exhibit  a  clear  1/A*  dependence. 


Fig.  3.  Schematic  representation  of  overall  scattering  geometry.  The  transmitter  is 
located  a  few  kilometers  from  the  lightning  activity.  Radio  signals  with  wave  vector  k 
are  incident  on  the  plasma  surface  at  an  angle  B  to  the  horizontal.  The  z  axis  is  taken  to 
along  the  axis  of  symmetry  of  the  plasma  channel,  which  is  approximately  perpendicular 
to  the  ground.  The  position  vector  s  extends  from  the  origin  of  the  coordinate  system  to 
the  surface  of  the  lightning-induced  plasma;  n  is  the  local  normal  vector  to  the  surface. 
h  is  the  vertical  extent  of  the  plasma,  which  is  estimated  to  be  on  the  order  of  kilometers. 
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Hg.  1.  Section  of  lightning-induced  plastiui.  showing  smooth  surface  irregularities 
with  characteristic  parameters  a  and  /  (figure  not  drawn  to  scale).  The  geometry  illus¬ 
trates  the  relation  k^R'  —  k^Ro  -  k2  •  s  used  after  Equation  (1). 


Fig.  5a.  Plot  showing  region  of  validity  of  1/A^  dependence  of  backscattered  power 
(Equation  (19)).  The  quantity  A/<t  must  lie  in  the  region  bounded  by  the  curves  47rcos0, 
■iTr{<r/l){cosd/ tand)  and  X/<r  =  10  for  the  approximations  made  in  the  analysis  to  be 
valid.  The  value  of  a/l  has  been  taken  to  be  0.1. 

Fig.  5b.  The  plot  of  Fig.  5a  is  repeated  for  «r//  =  0.05.  The  upper  bound  is  now 
given  by  A/cr  =  20,  which  is  above  the  point  of  intersection  of  the  lower  bounds. 
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Abstract 


The  problem  of  radio  backscatter  from  lightning-induced  pleismais  is  briefly  reviewed. 
The  predicted  wavelength  dependence  of  the  backscattered  power  is  compared  to  observa¬ 
tions  made  with  the  MIT  S-band,  C-band  cind  the  Millstone  UHF  radar.  The  analysis  is 
extended  to  plasma  disks  created  through  single  radar  beam  heating  of  the  atmosphere. 
The  artificial  plasmas  are  modeled  as  perfectly  conducting,  rough  plasma  pancakes,  with 
random  surface  inhomogeneities  described  by  either  by  Gaussian  or  power-law  type  spec¬ 
trum.  Only  the  limiting  case  where  the  angle  of  incidence  of  the  radio  waves  with  respect 
to  the  plasma  surface  is  large  is  considered.  In  radio  communication  and  space  surveillance 
applications,  the  reflectivity  of  the  plaisma  pancakes  may  be  independent  of  incident  wave¬ 
length,  if  the  surface  density  inhomogeneities  are  characterized  by  an  azimuthally  symmetric 
power  law  type  spectrum.  By  contrast,  a  wavelength  dependence  is  expected  in  diagnostic 
applications. 
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1  Introduction 


Radio  wave  reflection  from  the  ionosphere  has  long  been  used  to  extend  radio  communication 
paths  around  the  world.  The  daytime  F  peak  electron  densities  may  vary  between  10®-10®  cm“^, 
corresponding  to  plasma  frequencies  of  between  3-10  MHz.  The  ionosphere  may  therefore  reflect 
only  HF  and  lower  frequency  waves.  With  the  goal  of  extending  radio  communication  capabilities 
to  higher  frequencies,  radio  scientists  have  been  investigating  the  possibility  of  creating  high- 
density,  artificial  plasma  patches  in  the  atmosphere.  In  one  proposed  scheme,  two  very  powerful 
crossed  microwave  or  radio  wave  beams  are  employed  to  produce  thin  plasma  layers  through  air 
breaicdown  in  the  overlapping  beam  regions  [A.V.  Gurevich,  1980;  Kuo  et  al.,  1992].  In  other 
proposals,  a  single  radio  beam  is  used  to  “paint”  thin  plasma  patches  in  the  shape  of  the  beam 
cross-section  in  the  ionosphere  [Short  et  al.,  1991].  Transmitters  with  power  on  the  order  of 
gigawatts  may  produce  plasma  patches  dense  enough  to  reflect  radio  waves  with  frequencies  of 
up  to  1-2  GHz,  thus  extending  radio  communication  capabilities  by  two  orders  of  magnitude  in 
frequency  over  the  natural  ionosphere. 

Important  theoretical  considerations  arise  in  the  design  of  these  artificial  communication 
paths.  The  effects  of  surface  irregtdarities  on  the  reflectivity  of  the  plasma  patches  is  of  paxticulM 
interest,  since  they  may  impose  limitations  on  the  efficiency  of  radio  wave  reflection.  Scattering 
from  the  plasma  layers  created  in  the  crossed  beam  technique  is  essentially  a  Bragg  scattering 
process,  and  thus  unlikely  to  be  affected  significantly  by  the  small  surface  density  irregularities 
along  each  layer  [Kuo  et  al.,  1992].  By  contrast,  reflection  from  the  single  beam  produced 
“pancake”  plasmas  is  likely  to  be  significantly  affected  by  the  surface  roughness  of  the  plasma. 

Electric  storms  provide  a  natural  laboratory  for  testing  the  effects  of  surface  irregularities 


2 


on  the  radio  reflectivity  of  plasma  patches  in  the  atmosphere.  Lightning  is  a  series  of  short, 
complex,  pulsed  electric  discharges  that  create  highly  irregular,  dense  plasma  tubes.  Lightning 
backscatter  experiments  conducted  with  the  MIT  S-band,  C-band  and  the  Millstone  UHF  radar 
demonstrated  an  overall  inverse  square  dependence  of  the  radar  reflectivity  on  the  incident  wave¬ 
length  [Williams,  et  al,  1991].  This  behavior  may  be  successfully  explained  by  a  rough  surface 
scattering  model  of  the  lightning-induced  plasma  cylinders  [Dalkir,  et  al,]. 

This  paper  briefly  reviews  the  lightning  backscatter  experiments  and  the  rough  surface  scat¬ 
tering  model  that  may  account  for  the  inverse  square  wavelength  dependence  of  the  radio  wave 
reflectivity.  A  similar  model  of  rough  surface  scattering  is  then  applied  to  the  problem  of  radio  re¬ 
flection  from  single  radio  beam  induced,  pancake-shaped  plasmas.  Both  bistatic  and  monostatic 
radar  geometries  are  analyzed.  In  radio  communication  applications,  it  is  found  that  the  surface 
density  irregularities  of  the  artificially  induced  plasmas  may  introduce  a  wavelength  dependence 
into  the  radio  reflectivity,  if  the  spectrum  of  the  density  irregularities  is  Gaussian. 

This  paper  is  organized  as  follows:  Section  2  summarizes  the  MIT-MiUstone  lightning  backscat¬ 
ter  experiments;  the  theory  of  rough  surface  radio  wave  scattering  from  lightning-induced  plasma 
channels  is  outlined  in  Section  3.  Section  4  describes  the  application  of  the  anadysis  to  radztr 
scattering  from  single  beam  radio  wave  induced  plasma  disks.  Conclusions  follow  in  Section  5. 

2  Experimental  Results 

Lightning  backscatter  experiments  were  carried  out  with  the  S-band  (11  cm)  eind  C-band  (5.4 
cm)  radars  on  the  MIT  campus,  amd  the  steerablj  150  foot  UHF  (68  cm)  radar  dish  operated 
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by  the  MIT  Haystack  Observatory  at  Millstone  Hill  [Williams  et  al.,  1989;  Dalkir  et  al.,  1992]. 
The  radars  were  pulsed  for  1  to  3  //s  at  repetition  rates  of  between  200  to  900  Hz.  The  pulse 
echoes  were  recorded  over  a  wide  range,  from  a  few  kilometers  up  to  one  hundred,  over  different 
storm  days.  The  radar  beamwidth  was  1.4°  for  the  MIT  radars,  and  1.0°  for  the  Millstone  UHF 
radtir.  The  S-band  and  C-band  radars  were  operated  in  horizontally  polarized  mode,  while  the 
UHF  beam  was  circularly  polarized.  All  experimental  runs  were  conducted  with  the  radar  beams 
fixed. 

Figure  1  displays  the  reflectivity  of  the  plasma  at  the  three  different  wavelengths,  averaged 
over  numerous  individual  measurements.  The  mean  reflectivity  is  clearly  characterized  by  a  1/A^ 
dependence. 

3  Rough  Surface  Scattering  Model  of  Lightning  Backscat- 
ter 

This  section  outlines  the  theory  of  rough  surface  scatter,  as  applied  to  a  lightning-induced  plasma 
layer  in  the  atmosphere.  The  analysis  follows  the  method  introduced  by  M.  A.  Isakovich  (1952), 
and  later  extended  by  P.  Beckmann  and  A.  Spizzichino  (1963). 

Consider  the  scattering  of  plane  electromagnetic  waves  with  electric  field  where 

ki,  uj  are  the  wave  vector  and  angular  frequency,  respectively,  and  r  is  the  position  vector 
with  respect  to  the  origin  of  a  cylindrical  coordinate  system  with  basis  vectors  Cr,  e^,ex.  The 
scattering  plasma  surface  is  most  appropriately  modeled  as  a  long,  thin,  rough  cylinder,  with  its 
axis  of  symmetry  aligned  with  the  z-axis  (see  Figure  2).  The  vector  s  =  [5  -f  5s(z,  <^)]er  -I-  ze. 
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extends  from  the  origin  to  surface  of  the  plasma.  The  unperturbed  plasma  channel  has  radius 
S;  first-order,  reuidom  deviations  from  this  are  described  by  the  function  8s{z,  4>).  The  scattered 
wave  is  characterized  by  the  parameters  E2,k2  and  a;. 

In  this  simplified  model,  only  the  small-sc«de  density  irregularities  along  any  one  section  of 
an  approximately  vertical  lightning  channel  are  considered.  It  is  reasonable  to  assume  that  any 
density  perturbations  will  have  azimuthal  symmetry,  since  the  overdense  plasma  channels  are 
long  and  thin;  henceforward  we  assume  63(2,  (jy)  =  63(2). 

The  analysis  proceeds  by  use  of  the  Helmholtz  integral  which  describes  the  total  electric  field 
at  any  point  P  in  terms  of  the  electric  field  on  the  scattering  surface: 

where  /R',  R'  is  the  distance  from  P  to  the  scattering  point,  Ea  is  the  electric  field 

on  the  surface  A,  <uid  n  is  the  normal  vector  to  the  surface  at  the  scattering  point.  To  deal  with 
plane  waves,  we  consider  only  the  far  field,  i.e.  the  point  P  is  distant  from  the  scattering  point, 
such  that  k^R'  %  kiRo  -  k2  •  s  and  ^  «  exp(ife2i2o  -  *k2  •  s)/i2o  (see  Figure  3).  Furthermore, 
we  ignore  multiple  reflection  from  the  rough  surface,  so  that,  in  the  steady  state, 

Ea  =  {1  +  (2) 

=  i{l  -  R)E^ki  ■  (3) 

''On' -A 

where  R  is  the  reflection  coefficient  2ind  n  is  the  local  normal  to  the  surface  at  the  scattering 
point.  The  reflection  coefficient  is  in  general  a  function  of  the  eingle  of  incidence  and  polarization 
of  the  incident  waves,  and  the  electrical  properties  of  the  reflecting  material. 

Equation  (2)  is  a  good  approximation  to  the  surface  field  when  the  radius  of  curvature  of  the 
rough  surface  is  large  with  respect  to  the  wavelength  of  the  incident  wave,  i.e.  the  surface  has 
no  sharp  edges.  The  normal  derivative  of  the  surface  electric  field  is  obtained  by  considering  the 
reflected  and  incident  waves  with  respect  to  the  tangent  plane  to  the  surface  at  the  scattering 
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point.  Equation  (3)  is  valid  within  the  ray  approximation,  or  when  edge  (diffraction)  effects  can 
be  ignored.  This  is  the  case  for  L  »  X,  where  L  represents  the  characteristic  linear  size  of  the 
scattering  surface,  and  A  is  the  wavelength  of  the  incident  wave.  A  typical  lightning-induced 
plasma  cylinder  is  kilometers  in  length  and  may  be  meters  in  diameter,  and  thus  satisfies  the 
condition  L  >>  A  in  both  the  er  and  e,  directions  [Uman,  1969]. 

Substituting  Equations  (2)  and  (3)  into  the  Helmholtz  integral  (Equation  (1))  yields 

E{P)  =  1 1 [^(ki  -  k2)  •  n  -  (ki  -f  kz)  •  n]e'^^i-^2)  ‘^dA  (4) 

It  is  convenient  to  consider  the  normalized  electric  field  E  =  E/[iEie‘*^^/47riio].  Furthermore, 
kz  =  — ki,  since  monostatic,  diagnostic  radars  were  used  in  the  experiments.  Equation  (4)  can 
thus  be  rewritten  as: 


E  =  2rJ  Jk-  (5) 

where  k  =  kj;  R,  the  local  reflection  coefficient  at  each  scattering  point,  has  been  aproximated 

to  be  constant  over  the  entire  surface.  This  is  appropriate  for  the  lightning-induced  plasmas, 

which  are  almost  perfectly  conducting,  i.e.  (i2|  ~  1  [Williams  et  al.,  1989). 

The  surface  element  and  normal  vector  can  be  written  in  terms  of  the  density  perturbation 
5s{z)  eis: 

dA  =  Sd(t>dzy/l  +  {Ss'y  (6) 

n  =  {-Sa'eg  -f-  er)yH-  (5s')2 

where  6s'  =  d{Ss)/dz.  Substituting  these  terms  into  Equation  (5)  yields 

E  =  2R  r  Sd<t>  dz[-kM  +  (7) 

JO  J-hl2 

where  h  is  the  characteristic  height  of  the  scattering  cylinder. 
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The  mean  scattered  power  can  be  expressed  as  <  P  >=  (1/2)  <  EE*  >  /t/q  where  r/b  is 
the  intrinsic  impedance  of  the  propagation  medium;  rjo  =  1207r  ~  377  (17)  for  free  space.  From 
Equation  (7),  we  may  write 

<  EE*  >=  dzi  >  (g) 

k‘  J-h/2  J-h/2 

where  we  have  taken  i?*  =  1. 

To  proceed  further,  one  must  adopt  a  specific  model  for  the  roughness  parameter  Ss{z). 
Plasma  irrregularities  created  by  lightning  are  unlikely  to  have  a  periodic  or  similarly  determin¬ 
istic  structure;  rather,  noise-like,  random  perturbations  are  expected.  We  therefore  model  the 
extent  of  the  plasma  in  the  r  direction  as  a  random  variable  Ss(z),  which  is  completely  described 
by  a  Gaussian  joint  probability  density  function  /[5s(zi),  ^5(22)]  with  standard  deviation  <r,  and 
a  normalized  autocorrelation  function  C{\zi  —  22I)  with  characteristic  length  I  and  index  m: 


/[5s(2i),^s(22)] 


2'K(T^y 


exp(- 


6s{ziY  —  2C'C(.2i)^s(z2)  4-  8s{z2y 
20-2(1 (7*) 


Cdzi  -  Z2I)  =  exp(-lzi  -  Z2r//"') 


(10) 


Physically,  <t  and  I  characterize  the  height  and  size  of  the  irregularities  of  the  lightning- 
induced  plasma.  Any  roughness  profile  can  be  approximated  by  varying  the  parameters  of  the 
joint  probability  distribution  and  the  autocorrelation  function:  a  large  value  of  a  and  a  small  I 
imply  highly  localized  peaks  with  large  density  variations;  by  contrast,  small  o-  and  large  I  imply 
generally  low,  wide  irregularities  (see  Figure  3).  In  order  to  realistically  ignore  multiple  reflections 
in  the  model,  we  have  assiuned  the  absence  of  sharp  peaks  in  the  roughness  distribution  on  the 
surface  of  the  perfectly  conducting,  lightning-induced  plasma;  this  requirement  is  satisfied  in  the 
regime  I  »  X,  I  >  tr.  Previous  authors  have  concentrated  on  Gaussian  spectra  (i.e.  m  =  2  in 
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Equation  (10))  [Isakovidi,  1952;  Beckmann  and  Spizzichino,  1963];  studies  of  tr£insionospheric 
radio  wave  propagation  have  shown,  however,  that  power-law  type  spectra  provide  the  more 
realistic  model  for  plasma  irregularities  [Lee  et  al.,  1989,  cuid  references  therein].  We  therefore 
taJce  m  =  1  for  the  function  given  in  Equation  (10). 

The  definition  of  the  mean  in  Equation  (8)  can  then  be  rewritten  in  terms  of  the  probability 
density  functions  as: 

<  >=  r  dzi  r  dz2/(zi,Z2)e2*'^("^-"=')  s  (11) 

J —oo  •/— oo 

where  X2  is  the  two-dimensional  characteristic  function  of  the  joint  probability  distribution 
f[z\,Z2).  Similarly,  the  one-dimensioneil  characteristic  function  is  defined  as 

X{2K)  =<  >=  r  f{zy^’^‘dz  (12) 

J  ^OO 

where  f{z)  =  (V^<t)“^  exp(— z^/2<r^).  It  is  useful  to  first  calculate  the  variance,  d{E)=  < 
EE*  >  —  <  E  X  E  >*,  given  by 

d{E)  =  dz,  (iz2e'"''*<'‘-'=>[x2(2Av,  -2K)  -  x{2K)x\^K)]  (13) 

k‘  J-h/2  J-h/2 

For  a  Gaussian  joint  probability  distribution  function,  x(2fcr)  =  exp(— 2fc*<r*)  and  X2(2^,  — 2fcp)  = 
exp[— 4/!p<T^(l  —  C)];  thus  Equation  (13)  becomes 

d{E)  =  4^252  +  ^  dz,  (14) 

k‘  J-h/2  J-h/2 

The  integrand  vanishes  when  C  =  0,  or  when  |zi  —  22 1  is  large;  thus  we  may  appoximate  C  = 
exp(— |2i  —  Z2\) /I  ^  \  —  \zi  —  Z2\ll  and  replace  the  integration  limits  for  22  by  21  ±  /  >  22,  where 
I  is  the  maximum  range  of  |2i  —  22]  for  which  the  first  order  approximation  to  the  exponential 
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is  valid.  This  is  a  good  approximation  insofar  as  the  integrand  in  Equation  (14)  is  localized 
within  ±/;  it  can  be  seen  that  this  is  the  case  for  a  very  rough  surface,  i.e.:  >>  1,  which 

is  the  relevant  regime  for  the  lightning  backscatter  experiments  considered  here.  Integrating  the 
expression  in  Equation  (14)  and  neglecting  terms  on  the  order  of  yields  the  following  final 

result: 


<  P  >oc<  EE*  Z2hTT^S^ 


{kl  +  klf 


[i 


(15) 


kl  ^Akl^ 

The  term  {kl  +  kl)  =  k^  -  k^  indicates  that  <  P  >=  0  for  k  =  k^\  i.e.  the  backscattered  power 

is  zero  for  grazing  incidence  of  the  radio  waves  on  the  plasma  cylinders,  as  expected. 

The  result  may  be  applied  only  where  the  approximation  used  in  evaluating  the  integral  in 
Equation  (14)  is  justified;  this  is  the  case  for  Akl<r^  >>  1.  Multiple  reflections  from  the  scattering 
surface  may  be  ignored  only  insofar  as  the  density  irregularities  are  smooth;  this  condition  is 
met  when  I  »  \.  Furthermore,  if  >>  Ak^a^/P,  it  cein  be  shown  that 

hS^a^{kl  +  klY 


<  P  r 


2  7r2  I  kl 


oc  1/A^ 


(16) 


which  is  in  agreement  with  the  experimental  results  described  in  Section  2. 


4  Reflectivity  of  Single  Radio  Beam  Produced  Plasma 
Disks 

4.1  Diagnostics  with  a  Monostatic  Radar 

The  success  of  the  rough  sixrface  scattering  model  in  explaining  the  radio  reflectivity  of  lightning- 
induced  plasma  channels  suggests  that  it  may  be  usefully  employed  to  predict  the  behavior  of  a 


9 


radio  wave  induced  plasma  patch.  In  this  section  the  analysis  outlined  in  Section  3  is  applied  to 


the  pancake  or  disk-shaped  plasmas  produced  by  a  single  radio  beam. 

The  plasma  patch  is  modeled  as  a  disk  of  radius  S  and  height  h,  centered  at  the  origin  of  a 
cylindrical  coordinate  system.  Vectors  from  the  origin  to  the  the  bottom  surface  of  the  disk  axe 
given  by  s  =  rer  -h  [Sz{t^  <f))  —  /i/2]e55,  where  Sz{t,  <j))  is  the  random  variable  describing  the  height 
fluctuations  of  the  plasma,  characterized  by  the  Gaussian  joint  probability  distribution  given  in 
Equation  (9)  (See  Figure  4).  It  is  reasonable  to  assume  that  the  autocorrelation  function  of  Ss 
is  a  function  of  radius  only;  i.e.  C(|ri  —  r2|)  =  exp(— |ri  —  r2|'"//'”)  .  Thus  the  surface  density 
distribution  of  the  plasma  patches  can  be  taken  to  be  azimuthally  symmetric,  or  8z{r,  (f))  =  8z{r). 
Under  these  conditions,  one  may  write,  in  analogy  to  Equations  (6), 


n  =  (e,  -  +  (Sz'Y 

(17) 

dA  =  \Jl  +  {8z'Yrdrd<j> 

(18) 

where  8z'  =  d{8z)ldr.  Substituting  Equations  (17)  and  (18)  into  the  Helmholtz  integral  for  a 
monostatic  radar  (Equation  (5))  gives  the  following  expression  for  the  detected  electric  field: 

E  =  2^6-***''  r  d<t>  dr[k,  -  (ig) 

Jo  Jo 

Following  the  method  of  analysis  described  in  Section  3,  the  standard  deviation  of  the  normalized 
electric  field  is  found  to  be: 


d{E)  =  dr2  dn[Anr2  +  ^(rj  -  n)  +  (20) 

Jo  Jt^—1 

where  A  =  kl  +  2k^  +  ktlkl  =  k^/kl,  B  =  (Av/2i)[l and  G  =  k^/i^kt).  For  4A:f »  1, 
the  integrand  is  negligible  for  |ri  —  >  1.  Thus  we  may  make  the  transformation  p  =  —  r2 

and  estimate  the  integral  in  Equation  (20)  by 

d(E)  =  dr,  dp[Airl  +  p)  +  Bp  +  -  e*"''*"*]  (21) 

%  dr,  dp[A{r\)  +  Bp^ 
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Evaluating  the  integral  in  Equation  (21)  gives  the  following  expression  for  the  backscattered 


power: 


<  P  >oc<  EE"  167r^’ 


25W/Z  ,  ,  ,,,,,,  K^2kla^/l  , 


+ 


+  (5fc./8)[l  +  k^Jk% 


(22) 


^\2kt(T^IP  +  2k^  Akl  '  '  ^'{Aki<T^/P  +  klY' 

where  we  have  taJcen  iZ^  =  1,  and  where  a  power-law  type  of  irregularity  spectrum  [i.e.,  m  =  1 
in  the  autocorrelation  function  given  in  Equation  (10)]  has  been  assumed.  Terms  on  the  order 
of  exp{—4kla^)  have  been  neglected  in  Equation  (22).  We  consider  the  case  analogous  to  the 
lightning-induced  plasmas,  i.e.  the  surface  is  large  (S  »  A),  with  significant,  smooth  surface 
fluctuations  {1/<t  >  1,  Z  >>  A,  >>  1).  Under  these  conditions,  the  first  term  in  Equation 

(22)  is  dominant,  and  thus  the  backscattered  power  may  be  expressed  as: 

2ibVVZ 


<  P  >i 


(23) 


Hi  +  3*; 

This  result  is  basically  the  same  as  that  found  for  the  lightning  induced  plasmas,  given  by 
Equation  (15). 

If  a  Gaussian  autocorrelation  function  is  used  (i.e.  m  =  2  in  Equation  (10)),  the  backscattered 
power  is  given  by 


<  P  >«  v^(Z/2<T)()bVfc,^)(5V3/Z^)exp(-fc*ZV4fc*<r*) 


(24) 


While  convincing  evidence  shows  that  the  density  irregularities  of  the  lightning-induced  plasmas 
have  a  power-law  type  of  spectrum,  it  is  unclear  what  type  of  irregularity  spectrum  may  be 
produced  by  a  single  radio  beam  induced  plasma  patch.  In  the  czae  of  a  power-law  spectrum, 
the  radio  reflectivity  of  the  pancake  plasmas  is  expected  to  behave  analogously  to  the  lightning- 
induced  plasma;  by  contrast,  if  the  irregularities  are  more  properly  described  by  a  Gaussian 
spectrum,  the  radio  reflectivity  may  be  wavelength  dependent. 
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4.2  Radio  Communication/Space  Surveillance  with  a  Bistatic  Radar 

In  radio  communication  or  space  surveillance  applications  of  the  artificially  produced  plasmais, 
bistatic  radars  will  be  employed.  In  this  geometry,  the  incident  wave  vector  (ki)  and  the  detected, 
reflected  wave  vector  (k2)  are  described  by  ki  =  kr&r  +  K^x]  k2  =  k^Bx  —  kz^z-  The  k2  direction 
is  the  specular  direction  of  reflection. 

The  Helmholtz  integral  (Equation  (5))  thus  reduces  to 

E  =  2  y  jjRkz  +  (25) 

which  may  be  used  to  solve  for  the  power  reflected  in  the  k2  direction: 

d(E)  =  r'^'ir.Kr.r,  -ifc,/2(r,  +  -6-“!"“]  (26) 

•/O 

<P>oc  d{E)+  <E>  E>*=  (27) 

+  ^  -  2ie-^'’K5V3  +  klSHAk])]  +  S’*vV(4*f)] 

In  the  derivation  of  Equation  (27),  a  power-law  type  of  irregularity  spectrum  [i.e.  m  =  1  in 
Equation  (10)]  has  been  assumed.  Under  the  conditions  assumed  in  Section  3,  the  first  term 
dominates,  and  the  power  may  be  expressed  as: 

<p><x{i/6)^(s/,T)(i/<T)  m 

This  result  is  independent  of  both  the  wavelength  and  the  angle  of  incidence.  Thus  the  pancake 
plasmas  do  not  preferentially  scatter  any  one  wavelength  into  the  specular  direction,  if  the  plasma 
irregularities  have  power-law  spectra.  Equation  (28)  also  indicates  that  smoother  surfaces  (i.e., 
5//,  //<7  large)  scatter  more  power  into  the  specular  direction,  as  expected. 


12 


If  a  Gaussian  autocorrelation  function  is  assumed,  the  reflected  power  is  found  to  be 

<  P  >%  (l/6)V^fc,(5V^)//o-  (29) 

indicating  that  if  the  surface  density  irregulcirities  are  more  properly  described  by  a  Gaussian 
spectrum,  the  power  reflected  into  the  specular  direction  will  be  larger  for  longer  wavelengths 
at  any  given  angle  of  incidence.  As  mentioned  earlier,  it  is  unclear  whether  the  radio  wave 
produced  plasma  patch  has  a  power-law  or  Gaussian  type  of  plasma  irregularity  spectnim.  These 
theoretical  results  can  thus  serve  as  a  means  of  diagnosing  the  structure  of  the  artificial  plasma 
disk. 


5  Conclusions 

The  wavelength  dependence  of  the  backscattered  power  from  lightning-induced  plasmas  has 
been  studied  in  experiments  conducted  with  the  MIT  S-band,  C-band,  and  the  Millstone  UHF 
radars.  The  average  reflected  power  was  observed  to  decrease  with  the  wavelength  as  1/A^.  This 
behavior  can  be  successfully  explaiined  by  the  scattering  of  electromagnetic  waves  from  the  small- 
scaile  density  irregularities  of  the  lightning-induced,  perfectly  conducting  plasma  channels  in  the 
atmosphere.  These  density  irregularities  are  reeisonably  assumed  to  have  a  power-law  rather  thjm 
a  Gaussian  spectrum.  A  region  of  validity  for  the  predicted  1/A*  relation  is  believed  to  exist  in 
the  experiments,  which  depends  on  the  expected  size  of  the  density  irregularities  {o’ /I),  and  the 
angle  of  incidence  of  the  radio  waves  on  the  plasma  channels. 

Similar  results  axe  derived  in  the  case  of  monostatic  reflection  of  radio  waves  from  an  artifi¬ 
cially  produced  plasma  pancake.  In  bistatic  radio  communication/ space  surveillance  applications. 
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the  reflected  power  is  expected  to  have  no  wavelength  dependence,  if  the  plasma  irregularities 
have  a  power-law  spectrum.  By  contrast,  if  the  surface  density  irregularities  are  characterized 
by  a  Gaussian  spectrum,  the  reflected  power  is  inversely  proportional  to  the  incident  radio  wave¬ 
length.  Since  radio  communications  applications  of  artiflcial  plasma  patches  intend  primarily  to 
extend  ionospheric  reflection  to  higher  frequencies,  this  should  not  seriously  limit  the  efficiency 
of  the  proposed  schemes.  Thus  it  is  expected  that  radio  wave  induced  plasma  patches  may  be 
efficiently  employed  over  a  broad  range  of  frequencies  to  extend  radio  communication  paths  be¬ 
yond  their  present  capabilities 
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7  Figure  Captions 


Fig.  1.  Summ2iry  of  experimental  results.  The  average  volume  reflectivity  at  A  =  5.4,  11  and  68 
cm  is  plotted.  The  numbers  in  peirentheses  indicate  the  number  of  data  points  collected.  The 
wavelength  dependence  predicted  by  the  theory  of  long  thin  conductors,  and  underdense  plasma 
channels  is  also  indicated  on  the  figure.  Neither  existing  model  fits  the  data  of  the  present  study, 
which  exhibit  a  cleax  1/A^  dependence. 

Fig.  2.  Schematic  representation  of  overall  scattering  geometry.  The  transmitter  is  located  a 
few  kilometers  from  the  lightning  activity.  Radio  signals  with  wave  vector  k  are  incident  on  the 
plasma  surface  at  an  angle  6  to  the  horizontal.  The  z  axis  is  taken  to  along  the  axis  of  symmetry 
of  the  plasma  channel,  which  is  approximately  perpendicular  to  the  ground.  The  position  vector  s 
extends  from  the  origin  of  the  coordinate  system  to  the  surface  of  the  lightning-induced  plasma; 
n  is  the  local  normal  vector  to  the  surface,  h  is  the  vertical  extent  of  the  plasma,  which  is 
estimated  to  be  on  the  order  of  kilometers. 

Fig.  3.  Section  of  lightning-induced  plasma,  showing  smooth  surface  irregularities  with 
characteristic  parameters  a  and  I  (figvure  not  drawn  to  scale).  The  geometry  illustrates  the 
relation  k2R'  =  k2Ro  —  k2  •  s  used  after  Equation  (1). 

Fig.  4.  Geometry  of  radio  wave  scattering  from  pleisma  disks  of  radius  S  and  height  h.  The 
heater  array  "paints”  the  plasma  patch  through  air  breakdown  in  the  beam  cross-section.  The 
transmitted  wave  is  incident  on  the  plasma  at  the  origin  of  a  cylindrical  coordinate  system.  In 
the  diagnostic  geometry,  kbackacatter  =  — In  radio  communication  applications, 
the  reflected  wave  is  in  the  kgpecular  direction.  Note  that  the  scattering  surface  is  not  fiat,  but 
irregular  in  the  direction  (not  pictured). 
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